The effect of interactions on 2D fermionic symmetry-protected topological phases 

with Z 2 symmetry 



o 

Oh 
< 



i 

c3 



i 

c 
o 
o 



> 

On 

in 

o 
m 

> 



Zheng-Cheng Gu 1 and Michael Levin 2 

1 Institute for Quantum Information and Matter, 

California Institute of Technology, Pasadena, CA 91125, USA 

2 Condensed Matter Theory Center, Department of Physics, 

University of Maryland, College Park, Maryland 2074-2, USA 

(Dated: April 17, 2013) 

We study the effect of interactions on 2D fermionic symmetry-protected topological (SPT) phases 
using the recently proposed braiding statistics approach. We focus on a simple class of examples: 
superconductors with a Z2 Ising symmetry. Although these systems are classified by Z in the 
noninteracting limit, our results suggest that the classification collapses to Zs in the presence of 
interactions - consistent with previous work that analyzed the stability of the edge. Specifically, we 
show that there are at least 8 different types of Ising superconductors that cannot be adiabatically 
connected to one another, even in the presence of strong interactions. In addition, we prove that 
each of the 7 nontrivial superconductors have protected edge modes. 



Introduction. — Recently it has become apparent that 
generalizations of topological insulators pJE] known as 
"symmetry-protected topological (SPT) phases" [7HT5] 
can be realized in large classes of interacting boson and 
fermion systems. Loosely speaking, SPT phases are char- 
acterized by two properties. First, they support robust 
gapless boundary modes which are protected by certain 
symmetries. Second, SPT phases can be adiabatically 
connected to a "trivial state" (i.e., an atomic insulator 
or product state) if the relevant symmetries are bro- 
ken. While significant progress has been made in un- 
derstanding SPT phases in ID systems, [8) l9l H2HI6] less 
is known about the higher dimensional case. Several ap- 
proaches have been developed to understand these higher 
dimensional systems. One approach, which applies to 
bosonic SPT phases in general spatial dimension, is the 
cohomology classification scheme of Ref. [ITJ [T7] . An- 
other approach, which applies to bosonic or fermionic 2D 
SPT phases with chiral boson edge modes, is to study 
the edge theories of these systems using the if-matrix 
formalism [THH20]- 

In this paper, we discuss a third approach which was 
introduced in Ref. 5T! and applies to 2D SPT phases 
with unitary symmetry groups. The key idea behind this 
method is to study SPT phases by "gauging" their sym- 
metries - i.e. coupling them to an appropriate gauge 
field, thereby transforming their global symmetries into 
gauge symmetries. One can then probe the structure 
of the original SPT phases by constructing the excita- 
tions of the gauged systems and computing their quasi- 
particle braiding statistics. This approach has several 
nice features. First, it provides a simple way to distin- 
guish different SPT phases: if two gauged systems have 
different quasiparticle statistics then it is clear that the 
corresponding "ungauged" systems cannot be adiabati- 
cally connected without breaking the symmetry. Second, 
it gives insight into the stability of the edge: as shown 
in Ref. [51] the quasiparticle braiding statistics of the 



gauged system can be used to prove the existence of pro- 
tected edge modes. 

While Ref. HH focused on bosonic SPT phases, here 
we explore the fermionic case - a problem of particular 
interest because the classification of interacting fermionic 
SPT phases is not understood beyond ID (although an 
interesting attempt was made in Ref. l22|) . We focus on 
a simple class of examples: 2D superconductors with a 
Z2 Ising symmetry. It was previously conjectured[23 
|2"5] that while these systems are classified by an inte- 
ger invariant Z in the noninteracting limit, the classifica- 
tion collapses to Z§ when interactions are included. This 
claim was supported by an analysis of edge instabilities. 
Here we obtain further evidence supporting this conjec- 
ture. First, we show that there are at least 8 different 
types of Ising superconductors that cannot be adiabati- 
cally connected to one another, even in the presence of 
strong interactions. Second, we prove that each of the 7 
nontrivial superconductors have protected edge modes. 

Pseudospin notation. — We begin with some notation. 
Consider a general fermion system with an on-site, uni- 
tary Z2 symmetry S. Without loss of generality, we can 
assume that the Hamiltonian is built out of fermion op- 
erators that have a definite parity under S. [26] We will 
label the operators that are even under S with a pseu- 
dospin index f and operators that are odd under S with 
an index 4.. In this notation, the system is composed out 
of two species of fcrmions, c^ and cj,, where 



Sc+S 



c t; 



SciS 



-°i 



(1) 



In addition to the above Z<i symmetry, locality dictates 
that the system must also conserve fermion parity P/, 
defined by 

PfC^Pf = -^; PfC X Pj x = -c x (2) 

Putting these two constraints together, we can see that 
the pseudospin-f and 4- fermions are separately conserved 
modulo 2. 



The noninteracting limit. — We next review the clas- 
sification of noninteracting fermion SPT phases with Zi 
Ising symmetry. The key observation is that quadratic 
pseudospin mixing terms, e.g. clcj., are prohibited by 
the Z2 symmetry. Therefore the c^ and c^ fcrmions are 
completely decoupled in the non-interacting limit. Ap- 
plying the known integer classification of 2D topological 
superconductors [27l [28], it follows that the different free 
fermion phases are classified by a pair of integers (iA, iA). 



Here, 



y I 7/^1 



G Z corresponds to a phase where the 



pseudospin- ^ and pseudospin- 1 fermions form two decou- 
pled topological superconductors with iA and iA chiral 
Majorana edge modes, respectively. (The sign of iA,iA 
indicates the chirality of the edge mode - left or right 
moving) . 

In this paper, we only consider a subset of the above 
phases - namely those satisfying lA = —iA. The reason 
for this restriction is that our definition for SPT phases 
requires that they be adiabatically connected to a triv- 
ial band insulator if the symmetry is broken, and only 
phases with iA = —v^ obey this condition. Hence, ac- 
cording to our definition, the noninteracting SPT phases 
are classified by a single integer v = iA = — iA, 

The effect of interactions. — While the Z% symmetry 
requires that the pseudospin f and 4- fermions decou- 
ple from one another in the non-interacting limit, in- 
terspecies coupling is allowed once we add interactions 
into the system. (For example, the four fermion term 
cLcLcj^c^ is Zi symmetric, but mixes the two species). 
Thus, we might expect that the Z classification will col- 
lapse once we include interactions: i.e. it may be possible 
to adiabatically connect phases with different values of v. 
The question we will now investigate is: how many dis- 
tinct phases survive in the presence of interactions? 

For concreteness, we focus our analysis on a free 
fermion Hamiltonian with nearest neighbor(NN) hopping 
and pairing terms: 

W = zZzZ zZ ^ c La c ^;A + h.c. 

A=l <r (ij) 

V 

+ EEE %<r A ^ c L;A c i<r,A + h - c - ( 3 ) 

A=l <r (ij) 

Here i runs over lattice sites, while A = 1, ..., v describes 
different orbital states on each lattice site. We choose the 
the pairing term Ajjj-j-m so as to describe &p + ip (p — ip) 
superconductor, with an edge containing a single right 
(left) moving chiral Majorana mode. Our task is to deter- 
mine which H v describe distinct phases, and which can 
be adiabatically connected to one another. Our strategy 
for answering this question is to couple the pseudospin 
t and 4, fcrmions to two independent Z2 gauge fields, 



and therefore must represent distinct phases. 

The gauged model that we will analyze can be formally 
written as 

^gaugc = J2 £ £ Uj^T^Cj^ + h.C. 
A=l a (ij) 



EEE^U 



A 

ijcr jrr;X 



h.c. - H^ UX (4) 



A=l o (ij) 



where 



■ IJO 



is the Zi gauge field strength associ- 
ated with pseudospin a fermions and where H£ ux = 
E(yjw) T ij<r T jk<r T kla T L is a flux energy term that gives an 
energy cost to flux excitations of the gauge fields. (H^ ux 
is the analogue of the B 2 term in Maxwell electromag- 
netic dynamics). The Hamiltonian H gauge is defined in a 
Hilbert space consisting of gauge invariant states - that 
is, all states satisfying the constraint YIjpnnU) T ii<? = 
(— )^> ™i°-;*. This constraint can be thought of as a Z2 
analogue of Gauss' law, V • E = p. 

The next step is to compare the quasiparticle braiding 
statistics of the gauged system Eq. Q for different values 
of v. To this end, it is useful to first think about a simpler 
system with only one pseudospin component and v chiral 
edge modes. The quasiparticle braiding statistics of such 
a chiral superconductor were worked out by Kitaev in 
Ref. [551 That calculation showed that the quasiparticle 
braiding statistics of the superconductor depends on the 
number of chiral edge modes v, modulo 16. For example, 
if v is even, the Z 2 gauge fluxes (i.e. superconducting vor- 
tices) are Abelian anyons with an exchange phase factor 
ef w . If z^ is odd, the ^-fluxes are non- Abelian anyons 

(V 2 -1) 7T ■ 

with an exchange phase (— ) s e~s w when the two non- 
Abelian anyons are in the vacuum fusion channel. 

Now, let us consider the full system, which consists of a 
pseudospin j" component with v right moving edge modes 
and a pseudospin 4- component with v left moving edge 
modes. Naively, one might guess that the braiding statis- 
tics of this system also depends on v modulo 16, since it 
is made up of two independent chiral superconductors. 
However, this guess is incorrect: the braiding statistics 
of the "doubled" system only depends on v modulo 8. 
To see this, we need to show that the braiding statis- 
tics for v — 0,1,..., 7 are all different while the v = 
case is equivalent to the v — 8 case. One way to es- 
tablish the first statement is to compute the exchange 
phases of all the different types of Z% (or Z\) flux exci- 
tations. Here, a Z^ flux is defined to be any quasiparti- 
cle excitation that acquires a phase of —1 when braided 
around a pseudospin- \. fermion and acquires no phase 
when braided around a pseudospin-f fermion. Using the 
results of Ref. dHl it is easy to see that for even v there 



[Z\ x Zl) and then study the braiding statistics of the are 4 types of Z\ fluxes with exchange statistics ±e~> 



Z2 flux excitations in the gauged model. We will show 
that some values of v exhibit different braiding statistics 



while for odd v there are 2 types of Z^ fluxes with ex- 
change statistics ±e~^ ll/ . (In the latter case, we assume 



the fluxes are in the vacuum fusion channel) . In particu- 
lar, we see that the exchange statistics of the Z 2 fluxes are 
different for each of the eight possibilities v = 0, 1, ...,7. 
On the other hand, to see that v = and v — 8 have 
the same braiding statistics, we need to construct an 
explicit isomorphism between the quasiparticles in the 
two systems. To this end, we consult Ref. [35] and note 
that for both v = 0,8 the gauge theory Eq.Q has four 
quasiparticles 1, e CT , m CT , e a for each pseudospin direction, 
a = J [, 4- Including all possible composites of pseudospin 
t and 4- excitations, there are 4 • 4 = 16 quasiparticles 



all together. We can think of the e a as the constituent 
fermions while e a and m a are different types of Zf gauge 
fluxes which differ from one another by the addition of 
a fermion: e a — m a ■ e a . Using the results of Ref. (25] 
we can see that for both v = 0, 8, the three particles 
e a ,e a ,m a acquire a phase of —1 when braided around 
each other. The only difference is that e a and m a are 
bosons for the case v = while they are fermions for 
the case v = 8. With these properties in mind, one can 
easily see that the following map gives an isomorphism 
between the quasiparticles in the two systems: 



v = 1 e<f 
v = 8 1 6f£4. 



rwfei 






mi 



gtg+ 

m t mi 



myrni 






gtmj. m^-ei 



ete^ 



£te+ 



m t £4, gtrnj, 
m t mi 



r 



Here, the table is organized so that the first ten quasipar- 
ticles are all bosons while the other six are all fermions. 
We can see that the correspondence not only preserves 
braiding statistics and fusion rules, but also preserves 
the Z\ x Z^ gauge structure, mapping the 4- fermions 
(ei) of one system onto the corresponding fermions in 
the other system, and likewise mapping the Z\ fluxes 
(ej_, mj,, e^ei, £\rni) of one system onto the Z^ fluxes of 
the other system (and similarly for f). 

Two conclusions follow from the above analysis. First, 
we conclude that the Hamiltonians H v with v = 
0, 1, • • • ,7 cannot be adiabatically connected to one an- 
other without breaking the Zi symmetry. Indeed, if there 
existed a gapped path connecting these Hamiltonians, 
then there would have to be a corresponding path con- 
necting the gauged systems H sausc - an impossibility, 
since we have seen that they have different quasiparticle 
braiding statistics. The second conclusion is that it is 
at least plausible that H° and H 8 can be adiabatically 
connected to one another in the presence of interactions, 
since the corresponding Z 2 x Z\ gauge theories share the 
same statistics and gauge structure. 

The instability ofv = 8 edge. — In this section, we give 
additional evidence that the v = 8 system is a trivial 
SPT phase: we show that the v = 8 edge can be gapped 
out by appropriate interactions, without breaking the Z 2 
symmetry (explicitly or spontaneously). We note that a 
similar result was obtained previously in Refs. E5H251 

Our approach is based on bosonization. We note 
that the edge of the v — 8 free fermion Hamiltonian 
Eq.([3]) contains 8 pseudospin-f Majorana modes and 8 
pseudospin-4 Majorana modes moving in opposite direc- 
tions. Pairing up the Majorana modes to form com- 
plex fermions, we can equivalently describe the edge us- 
ing 4 pseudospin-^ and 4 pseudospin-4 complex fermions. 
We then bosonize these fermions, using 4 boson modes 
<&!,.. ,$4 for the pseudospin-^ fermions, and 4 boson 



modes $5, ..., $g for the pseudospin-4 fermions. The edge 
is then described by the chiral boson Lagrangian 

Cedge = ^(Kjjd^jd^j - V U d x *id x *j) (5) 

where K = diag(l, 1, 1, 1, — 1, — 1, — 1, — 1), and Vjj is 

the velocity matrix. Here we use a normalization con- 
vention where the fermion creation operators are of the 
form e l * fc , k = 1,...,8. In this language, the symmetry 
transformation is given by S~ 1 &S = $ + irK ~ l \ where 

x 1 



(0,0,0,0,1,1,1,1) 



We now construct interaction terms that gap out the 
edge without breaking the Zi symmetry (either explicitly 
or spontaneously). We consider backscattering terms of 
the form U(A) = U(x)cos(A T K& — a(x)). In order for 
[/(A) to be invariant under 5 1 , we require that 



A 1 X = (mod 2) 



(6) 



In order to gap out the edge, we need to add 4 backscat- 
tering terms J2i U(Ai): each term can gap out a pair of 
counter-propagating edge modes. Such terms can gap 
out the edge as long as the {A^} vectors satisfy [32] 



aTka, = o 



(7) 



for all i,j. This "null- vector" condition guarantees that 
we can make a suitable change of variables mapping L e d ge 
onto a system of 4 decoupled Luttinger liquids with 4 
backscattering terms. It is then easy to see that the 
backscattering terms will gap out the corresponding Lut- 
tinger liquids (at least for large U [5TJV 

We now claim that the following {A,} will do the job: 

Aj = (1,-1,0,0,1,-1,0,0) 



Al = (1,0,-1,0,1,0,-1,0) 
Aj = (1,0, 0,-1, 1,0, 0,-1) 
Aj = (1,0,1,0,0,-1,0,-1) 



(8) 



U a U b 
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(b) 




W f 




W f 



Wa 



(c) 



(d) 



FIG. 1. (a) We consider a thought experiment in which we 
create two Z2 fluxes in the bulk and then move them along a 
path p to points a, b at the edge, (b) We argue that the two 
fluxes can be annihilated at the boundary by applying local 
operators U a , Ub- (c) We define W/3 to be an operator which 
describes a process in which the fluxes are created in the bulk, 
brought to the edge, and then annihilated, (d) To obtain a 
contradiction, we consider two paths P, 7 that intersect one 
another, and we investigate the commutation algebra of the 
corresponding operators W,g,W 7 . 

Indeed, it is easy to check that these {A^} obey the null 
vector criterion d7b, as well as the symmetry condition 
(pi). To complete the argument, we need to check that 
the perturbation corresponding to {A^} does not sponta- 
neously break the Zi symmetry. However, as explained 
in Ref. [18], we can rule out the possibility of sponta- 

neous symmetry breaking if the 4x4 minors of the 

8x4 matrix with columns Ai , . . . , A4 have no common fac- 
tor. This property of Ai, ..., A4 can be verified by direct 
calculation. 

Protected edge states for v 7^ mod 8 . — On the other 
hand, we now show that the edge of H v is protected if 
v 7^ mod 8. To state our result more precisely, let us 
consider a disk geometry and a Hamiltonian of the form 
H = Hbulk + -Hodge, where iJ bu ik = H v , and H cdgc is 
an arbitrary interacting Hamiltonian acting on fermions 
near the edge. In this setup, what we will show is that 
the ground state |0) cannot be both Zi symmetric and 
"short-range entangled." [33] We believe that this result 
rules out the possibility of a Z^ symmetric, gapped edge, 
and in this sense proves that the gapless edge excitations 
are protected. 

As in Ref. [3T] , our argument is a proof by contradic- 
tion: we assume that |0) is short-range entangled and Z2 
symmetric and we show that these assumptions lead to a 
contradiction. The first step is to couple the pseudospin-^ 
and pseudospin-j, fermions to two independent Z2 gauge 



fields as in Ml). We then imagine creating a pair of Z\ 
(or Z\~) fluxes in the bulk. After creating the Z% fluxes, 
we separate them and move them along some path /3 to 
points a, b at the boundary (Fig. [Ik). Formally, this pro- 
cess can be implemented by applying a unitary (string- 
like) operator Wp to |0). 



Next, we claim that the Z\ fluxes can be annihilated 
at the boundary if we apply appropriate local operators. 
That is, there exist local operators U a ,Uh, acting near 
points a,b such that U a U b Wp\0) = |0) (Fig. ^p). Es- 
tablishing this claim is the hardest step in the argument, 
and here we merely outline its proof. [33J The basic point 
is that when we bring the Z2 fluxes to the boundary, 
we effectively create two Z^ domain walls at a and b. 
Given that the ground state is Z\ symmetric, these do- 
main walls are local excitations: they only affect expecta- 
tion values in the neighborhood of a and b. It then follows 
that these domain walls can be annihilated by local op- 
erators since local excitations of a short-range entangled 
state can always be annihilated locally. 

In the third step, we consider a creation and annihi- 
lation process in which two Z\ fluxes are created in the 
bulk, moved to the boundary and then annihilated. Let 
W^3 be a unitary operator describing this process (Fig. 
fit). (Formally, Wp = UJJbWp). Now, consider a sec- 
ond path 7 with the geometry shown in Fig. [TH and 



define 



in the same way. By construction, we have 



W/}|0) =W T |0) = |0). Hence 



,|0) =W 7 W /9 |0) = |0) 



(9) 



In the final step, we show that |9]) leads to a contradic- 
tion if v 7^ (mod 8) . It is useful to consider separately 
the case where v is even and v is odd. First, suppose v is 
even. In this case, the Z\ fluxes are abelian anyons and 
it follows from a general analysis of abelian quasiparticle 
statistics (see e.g. Refs. I!?l"l and I34"|) that 



W/3W 7 |0) 



Jlid 



W 7 W^|0) 



(10) 



where e l8 is the exchange phase of the Z\ fluxes. Ac- 
cording to the braiding statistics calculation outlined 
above, the four types of Z\ fluxes have exchange statis- 
tics = ±Tf. Hence if v ^ (mod 8) then e 2W ^ 1 for 
any of the four types of fluxes and equations ([9]), ( 10 ) are 
in contradiction. 

Now suppose v is odd. In this case, the Z^ fluxes are 
non-abelian anyons, so the above braiding statistics anal- 
ysis is more complicated. However, we can avoid these 
complications using an alternative argument. We note 
that if v is odd then each Z\ flux carries an unpaired 
Majorana mode. Thus, the state W^|0) has unpaired Ma- 
jorana modes localized near points a and b. But then it 
is clearly impossible for U a U~bWp\0) = |0) since unpaired 
Majorana modes cannot be destroyed by any local opera- 
tion. Once again, we encounter a contradiction, implying 
that our assumption is false and |0) cannot be both Z2 
symmetric and short-range entangled. 

Conclusion. — In this paper we have studied SPT 
phases in interacting fermion systems using a braiding 
statistics approach. As a simple example, we considered 
superconductors with a Zi (Ising) symmetry. Although 
in the noninteracting case these Ising superconductors 



are classified by an integer invariant v £ Z, we give evi- 
dence that the classification collapses to Z§ in the pres- 
ence of interactions. We also give a general argument 
proving that the edge excitations are protected when 
v 7^ (mod 8) and unprotected when v — (mod 8). 
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